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1 Abstract 

We consider L 2 minimizing geodesies along the group of volume preserving 
maps SDif f(D) of a given 3-dimensional domain D. The corresponding 
curves describe the motion of an ideal incompressible fluid inside D and 
are (formally) solutions of the Euler equations. It is known that there is a 
unique possible pressure gradient for these curves whenever their end points 
are fixed. In addition, this pressure field has a limited but unconditional 
(internal) regularity. The present paper completes these results by show- 
ing: 1) the uniqueness property can be viewed as an infinite dimensional 
phenomenon (related to the possibility of relaxing the corresponding mini- 
mization problem by convex optimization) , which is false for finite dimen- 
sional configuration spaces such as 0(3) for the motion of rigid bodies; 2) 
the unconditional partial regularity is necessarily limited. 
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2 Introduction 

There are very few problems in mathematical fluid mechanics for which global 
existence and uniqueness results are known without any restriction on the 
data. (For instance, the Leray global existence result for 3D Navier-Stokes 
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equations lacks uniqueness, the Yudovich theorem for 2D Euler equations re- 
quires bounded vorticity, the Glimm-Bressan theory for ID gas dynamics is 
only valid for small initial data, etc.). One of them is the problem of min- 
imizing geodesies ('shortest paths') on the group of 3D volume preserving 
diffeomorphisms with L 2 metric, which, following Arnold's geometric inter- 
pretation |AK] of the Euler equations, is a secondary way to find solutions 
of these equations without solving an initial value problem. Let us be more 
specific. 

Given a smooth enough bounded domain D in R 3 , we denote by SDiff(D) 
the group of all volume and orientation preserving diffeomorphisms of D. 
This group is a natural configuration space for the motion of an incom- 
pressible fluid moving inside D which is mathematically described by a time 
dependent curve t — > gt G SDif f(D). The minimizing geodesic (or shortest 
path) problem can be defined as follows: given two such diffeomorphisms 
go, g\, find a curve t G [0, 1] — >• gt G SDiff(D) that achieves the geodesic 
distance between g and g\. 

Ato,»)=irf ) / J £ 1 ||f|li,«« ID 

where || • | denotes the norm in L 2 (D, R 3 ). The formal optimality equation 
reads 

^oft-i + V ft = 0, (2) 

where p t is a time dependent scalar field defined on D (called the 'pressure 
field') which balances the incompressibility constraint as a Lagrange multi- 
plier. Equation (j2J) is precisely the Euler equation introduced in 1755 |Euj 
to describe the motion of an inviscid incompressible fluid moving inside D 
without any external force. Let us empasize that the minimizing geodesic 
problem is different from the more conventional Cauchy problem, for which 
the initial 'velocity' %r is prescribed together with g (which is traditionnaly 
normalized to be the identity map) and, of course, there is no prescribed 
endpoint g\ at time t = 1. The minimizing geodesic (or shortest path) prob- 
lem has been solved at the local level (together with the Cauchy problem) by 
Ebin and Marsden |EMj : if g\ o g^ 1 — J is sufficiently small in Sobolev norm 
H s , s > 5/2, there is a unique minimizing geodesic. In sharp contrast, a 
striking result of Shnirelman |Shl] shows that the minimizing geodesic prob- 
lem may have no solution at all in the large. In the same paper, it is proven 



2 



The minimizing geodesic problem in fluid mechanics 



that, as a metric space, the completion of SDiff(D) for the geodesic dis- 
tance (OQ) is the semi-group VPM(D) of all volume-preserving maps g of D 
in the measure theoretic sense: g is Borel and 



f(g(x))dx = / f(x)dx, 

D JD 

holds true for all continuous functions /. Notice that, in particular, any 
map g(x) = {h{x\), x%, X3) , for which h is a Lebesgue measure-preserving 
map of the unit interval [0, 1], belongs to the completion of SDiff(D) when 
D = [0, l] 3 ! However, even in the completed configuration space VPM(D), 
the minimizing geodesic problem may have no solution, as shown by the 
author in [Brl] with the example: 

D = [0,1] 3 , g (x), gt(x) = (h(x 1 ),x 2 ,x 3 ), h{x 1 ) = 1 - xi. 

More positively, it was proven in [Br2l IBr3] that there is a unique pressure 
field attached to the minimizing geodesic problem. With the help of a crucial 
density result obtained by Shnirelman [Sh2j . this result can be stated in the 
following way: 

Theorem 2.1 Let g and g\ given in VPM(D) where D = [0, l] 3 . We say 
that g\ G SDiff(D) is an approximate minimizing geodesic between g and 
9i if 

£\\^jjr\\hdt^ A(g , gi ) 2 , \\g e -g Q \\l>^0, M - gi\\h -» 0, (3) 

as e l 0. Then, there is a unique pressure gradient field Vpt such that, for 
all approximate minimizing geodesic, 

holds true in the sense of distributions in the interior of [0, 1] x D. 

The main idea of the proof is that the minimizing geodesic problem, which 
apparently is a minimization problem lacking both compactness and convex- 
ity, can be relaxed as a CONVEX minimization problem, in an appropriate 
generalized framework, without relaxation gap, thanks to Shnirelman's den- 
sity result. This is strongly related to the fact that the completed configura- 
tion space VPM(D) is a dense subspace (as shown in [Net IBGj ) of a weakly 
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compact convex set, namely DS(D), the set of all doubly stochastic mea- 
sures on D x D (i.e. all nonnegative measures on D x D having the Lebesgue 
measure as projection on both copies of D), with respect to the weak conver- 
gence of measures, through the embedding g 6 VPM(D) — > \i g 6 DS(D), 
where 



The uniqueness of p follows, quite easily, from a duality argument due to 
this 'hidden' convex structure, as shown in |Br2j . Nevertheless, geometri- 
cally speaking, the uniqueness of p is quite surprising, since, between two 
given points, minimizing geodesies are not necessarily unique (as can be eas- 
ily checked) whenever they exist and there are no a priori reasons that the 
corresponding acceleration fields — Vp should be identical. It is unlikely that 
such a property could be proven using classical differential geometric tools. 
Theorem 12.11 can be completed in various manners, as in |Br3] (existence of 
generalized minimizing geodesies solving a generalized version of the Euler 
equations, partial regularity of p, etc.) or, recently, by Ambrosio and Figalli 
in [AF14 IAF2] . Bernot, Figalli, Santambrogio [BFSJ. In particular, Vp is 
shown to be a locally bounded measure in the interior of [0, 1] x D in [Br3j. 
The time integrability of p has been improved since by Ambrosio and Figalli 
|AF2j : p belongs to the space Lf oc (}0, 1[, BVi oc (D°)). At this point, two nat- 
ural questions can be asked: 

Ql: Is the uniqueness of the pressure a specific property of minimizing 
geodesies on the infinite dimensional group of volume preserving diffeomor- 
phisms, which is not true for similar finite dimensional configuration spaces, 
such as the special orthogonal group SO (3) for the motion of a rigid body? 
Q2: Is the regularity of the pressure field only partial? 
In both cases, the present paper provides a positive answer. 

3 About the uniqueness of the pressure field 

In this section, we give some evidence that the uniqueness of the pressure 
field for the minimizing geodesic problem on SDif f(D) (or its completion 
VPM(D)) is a genuine infinite-dimensional phenomenon. For this purpose, 
we consider the finite dimensional situation where the special orthogonal 
group 5*0(3) substitutes for SDiff(D) and rigid motions of solid bodies 
substitute for incompressible inviscid fluid motions. (See |Arj.) As a matter 
of fact, as shown below, rigid motions can be interpreted as particular solu- 
tions of the Euler equations, in the special case when the fluid domain D is 
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an ellipsoid D = KBi, where B\ is the unit ball in R 3 and K a symmetric 
positive matrix. To see that, we first notice that a volume preserving map 
g : D = KBi -)-Disa linear map if and only if: 

g(x) = KUK~ l x, xeD = KBi, (4) 

for some matrix U in the orthogonal group 0(3). Indeed, by definition of 
D = KBi, x — » K~ x g{Kx) must be a linear map preserving the unit ball 
£>!, i.e. an element of 0(3). The condition U G 0(3) can be expressed in a 
variational way by requiring 

Tr(U*MU) = Tr(M) (5) 

for all 3 x 3 real symmetric matrices M, where Tr and * respectively denote 
the trace and the transposition operator on 3 x 3 real matrices. 
Let us now look at geodesies (g t ) on SDiff(D) with the additional constraint 
that gt must be linear at each time t ('rigid motions'), i.e. of form fl4]): 
gt = KUtK~ l for some curve Ut valued in 0(3). We have 

[ \^fdx= [ \K^K-^dx 
Jd dt Jd dt 

Jb 1 dt dt dt 

where r , r\ > are normalization factors. Thus, encoding the condition 
U G 0(3), the geodesies we are looking for correspond to saddle points 
(U t ,M t ; t G [0, 1]) of the Lagrangian 

j o {Tr(^-K 2 ^) - Tr(U t *M t Ut) + tr(M t )}dt. (6) 

Here M t is constrained to be symmetric and there is no more constraint on 
Ut- The optimality conditions are straightforward: 

K^ + M t U t = 0, (7) 

with U t orthogonal and M t symmetric. At this point, we have exactly recov- 
ered the usual equations for motions of a rigid body in classical Mechanics 
(see [Arj). They describe geodesic curves on 5*0(3) for the metric generated 
by the 'inertia' matrix K 2 of the body. 



5 



The minimizing geodesic problem in fluid mechanics 

The point of our discussion is that rigid motions not only can be embedded 
in the framework of fluid Mechanics but, more strikingly, are just special 
solutions of the Euler equations. Indeed, from ([7]), we recover a solution of 
the Euler equations just by setting 

g t (x) = KU t K- l x, p t (x) = -K- x M % K' x x -x, xeD = KB 1 . (8) 

We conclude that geodesies on 5*0(3) are just special geodesies on SDif f(D) 
for a special choice of the domain D (namely the ellipsoid D = KB\ associ- 
ated to the inertia matrix K 2 ). 

Now, we can discuss the issue of minimizing, not just plain, geodesies. The 
key point is that, for a solution of the rigid body motion (J7J) and the associated 
solution OH]) of the Euler equations, which corresponds to a geodesic curve 
on both 50 (3) and SDiff(D), it is NOT equivalent to be a minimizing 
geodesic on 50(3) and SDif f(D). This makes sense, since rigid motions are 
more restrained than fluid motion. So it is conceivable that Ut is a minimizing 
geodesic on 5*0(3), meanwhile the corresponding g t = KUtK x is a plain, 
non minimizing, geodesic on SDif f(D). This suggests that the uniqueness 
of the pressure field for fluid motions has no equivalent property in the case of 
rigid motions. As a matter of fact, this follows from the following (stronger) 
result: 

Theorem 3.1 Let (Ut,M t ), (Vt,N t ) two distinct solutions of the rigid mo- 
tion equations for t G [0, 1] such that Uq = Vq = I (where I is the identity 
matrix) and U\ = V\. Assume that M t = N t for all t e [0, 1]. Then Ut and 
Vt are exceptional in the sense that they must be rigid rotations with constant 
angular speed around one of the inertia axis of the body (i.e. an eigenvector 
of the inertia matrix K 2 ). 

In other words, generically, two distinct geodesies meeting at two different 
points must have different accelerations. This is, a fortiori, also true for 
minimizing geodesies. The group 50(3) is a smooth closed bounded set in 
the finite dimensional Hilbert space of real 3x3 matrices and has a finite 
geodesic diameter. Thus, every geodesic curve can be minimizing only on 
finite time intervals, and, therefore, they are plenty of distinct minimizing 
geodesies connecting a same pair of points. 
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Proof of theorem 13.11 

Let us consider two distinct minimizing geodesies on SO(3), (U t , M t ), (V t , N t ), 
with same endpoints U = I, U\ at t = and t = 1. It is classical to introduce 
two fields of skew symmetric real matrices, B t and C t , such that 

^ = B t U t , ^ = C t U t 
dt ' dt 

(this is like introducing Eulerian coordinates in fluid Mechanics). Then, the 
motion equation ([7]) can be written 

^ + B\ = -K~ 2 M t , (9) 

and, similarly 

^ + C 2 = -K~*N, 

Let us assume that M t and N t coincide, while B t and Ct are distinct. Since 
Bl and C 2 are symmetric meanwhile B t and Ct are skew symmetric, it follows 
that Ct = B t + L for some constant skew symmetric matrix L different from 
zero, meanwhile C\ = _B t 2 . This implies B t L + LB t + L 2 = 0. and, therefore, 

dB t r dB t 
—±L + L—± =0. 
dt dt 

Up to a change of orthonormal frame, we may assume, without loss of gen- 
erality, Lx = (x 2 , — Xi, Q)f3 for some constant (3 ^ 0. We deduce, by di- 
rect calculation, that ^ = 0. So B t = B , C t = B + L. Using again 
B t L + LB t + L 2 = 0. we get no other solution than B x = (—x 2 ,Xi,0)/3/2 
and C = —Bq. This means that U t = exp(_B t) and V t = exp(—B Q t) are just 
rotations at constant angular speed +(3/2 and —(3/2, along the axis (0, 0, 1). 
Notice that, since the time interval for which the geodesies are minimizing 
has been fixed to be [0, 1], the only possibility is (3 2 = An 2 . Going back to 
(El), we further deduce M t x = M x = —K 2 BqX = K(xi,X2,0)n 2 where M t 
is supposed to be symmetric. This implies that the axis of rotation (0, 0, 1) 
must be an eigenvector for the inertia matrix K 2 (i.e. an inertia axis for 
the rigid body), which is clearly an exceptional situation, as soon as K is 
a generic symmetric positive matrix (which corresponds to a generic rigid 
ellipsoid). 
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4 Limited regularity of the pressure field 

In this section, we provide an explicit, self-similar, solution of the minimizing 
geodesies with limited regularity. 

Theorem 4.1 Let L > 0, D = [-L, L] x [0, If, L>1. 

Then, for < t < 1, x E D, 

g t (x) = t 2/3 (2 v /x 1 t- 2 / 3 - 1, z 2 , x 3 ), < x x < t 2/3 , (10) 

g t (x) = t 2/3 (l - 2^- Xl t- 2 / 3 ,x 2 ,x 3 ), -t 2/3 <X!<0, (11) 
g t (x)=x, |xi|>t 2/3 , (12) 

define a (generalized) minimizing geodesic, with a pressure field of limited 
regularity 

p t (x)=p t (x 1 ) = -^{t^-x\) + . (13) 

Remarks 

i) The family (g t , < t < 1) is not valued in SDif f(D) but in its completion 
VPM(D). So, our example does not prevent a better regularity of the pres- 
sure field in the case of smooth data g , g±, valued in SDif f(D). However, 
it does rule out unlimited internal regularity, in the style of classical elliptic 
PDE theory, independently on the boundary data. This example also shows, 
in our opinion, that the regularity to be expected for the pressure field is 
semi- concavity, or, at least, measure- valued second order space derivatives. 
Let us recall that, so far, we only know that Vp is a locally bounded measure 
[Br5] (or more precisely p belongs to Lf oc (}0, 1[, BV loc (D°)) [5F2] ). So there 
should be one order of differentiability in space to be gained in the future. 

ii) This solution can be interpreted as a 'hydrostatic vortex sheet' as ex- 
plained in section and, under that form, coincides with a self-similar 're- 
laxed solution' of the Euler equations already introduced by Duchon and 
Robert in jDRj . 

Proof 

The proof is based on two statements: 

i) For each t, gt is a volume preserving map (which is not obvious at first 
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glance). 

ii) For almost every fixed x G D, £ t = gt(x) minimizes 

//{i|f I 2 - ta)(6)K (14) 

among all curve £ t valued in D such that £o = 9o( x ), £1 = 9i( x )- These two 
conditions guarantee that, indeed, g% is a minimizing geodesic, following a 
standard argument (see |Br3l IAFT] . for instance). 

Let us first check the second statement. If \x-±\ > 1, gt{x) = x and the 
statement is trivial. If \xi\ < 1, and x\ 7^ 0, we see that £ t = gt(x) is 
continuously differentiable in t G [0,1], with piecewise continuous second 
order derivative (with a jump at t = |xi| 3//2 ), and satisfies 

^ = -(V Pt )(6) (15) 

for almost every t, which guarantees that £ is already a critical point of the 
action defined by (TH1) . Let us now prove that £ is also a global minimizer 
of the action (TH1) . For this purpose, we compute the second variation SV 
of action f[T4"|) for a perturbation £ t + with Co = Ci = 0, and want to show 
that SV > 0. We find 

^ = i{/i§i 8 -(9 a p t )((e*)i)(c«)?K 



By definition (fT3l) of p 4 
2 

9I2 



d2 Pt{s) = £;[\{\8\ < t 2/3 } - t 2/3 5( S - t 2/3 ) - t^5(s + t 2/3 )] 



Thus 



2 Jo dt ] 9t 2 7|(^)i|<* 2 / 3 
Let us recall the classical Hardy inequality (in one space dimension): 

for all smooth real function r\ t such that 770 = 0. We deduce that SV > 
and conclude that £, indeed, is a minimizer of (1141) as its end points are fixed 
at t — and £ = 1. 
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Finally, let us prove that g% is volume preserving. Due to the self-similarity 
of gt, it is enough to check, in the case t — 1, that, for every real continuous 
function /, 

/ f(l - 2V^x)dx + f 1 f(2y/x - l)dx = f 1 f(x)dx 

J-l JO J-l 

which also follows from elementary calculations. 

Finally, let us mention that the solution discussed in this section has a fluid 
mechanic interpretation and can be derived from the hydrostatic limit of the 
Euler equations. See all details in section |5j 

5 Appendix on the hydrostatic limit of the Euler equa- 
tions 

As explained in |Br6] . the minimizing geodesic problem is strongly linked 
to the hydrostatic limit of the Euler equations, which reads, on the domain 
D=[-L,L] x [0,1] 2 , 

D tVl + d lP = 0, D t v 2 + d 2 p = 0, D t = d t + v-V, (17) 

3 p = 0, V-v = 0, v//dD. (18) 

These equations are formally obtained by ignoring the vertical acceleration 
term D t v 3 in the classical Euler equations (see (LTJ, IBr4| IGr[ IBr5] for some 
rigorous results). Notice that, given a sufficiently smooth solution (v,p) 
of these hydrostatic equations, we may introduce the corresponding flow 
X(t,x), defined by 

d t X(t,x) = v(t,X(t,x)), X(t = 0,x)=x, (19) 

which provides a time dependent family of volume-preserving maps X(t, •) 
of D, since v is divergence- free and parallel to the boundary. We are going 
to construct an explicit solution to these equations and, as an output, the 
solution flMTITT2lfT3ri used in Theorem EH 

Construction of an explicit solution 

We first define a divergence-free velocity field (with trivial second component) 
by setting: t> 2 = 0, v i = d 3 ip, v 3 = —diip, where ip is the 'stream-function' de- 
fined by three different formulae in the domain D, depending on the location. 
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We first set: 

1>(t t x 1 ,Xt) = rl/ * Xs jf~ 1 \ 0<x 3 <^, |e|<l, (20) 
where we use the rescaled coordinate £ = x\t~ 2 ^. Next: 

m^x 3 ) = t ' 1/3{X3 ~ 3 m + 1 \ i±l<xs<l, |£|<1, (21) 
and, finally: 

il>(t, Xl ,x 3 )=0, |£|>1. (22) 

Notice that the stream-function is continuous at the interfaces £ = 1, £ = 
— 1 and x 3 = (1 + £)/2 and we can easily recover the velocity field t> by 
differentiating -0: 

«i = <9 3 ^ = , v 2 = 0, v 3 = -dxij) = — (23) 

whenever < x 3 < |£| < 1, 

"l = g , ^2 = 0, w 3 = (24) 

whenever ^ < £3 < 1, |£| < 1 , and t> = whenever |£| > 1. 
This velocity field is piecewise smooth, with a strong singularity at t — and 
also at the interfaces £ = 1, £ = —1, x 3 = (1 + £)/2 for each t > 0. The 
interface x 3 = (1 + £)/2 can be interpreted as a vortex sheet initially located 
verically above x\ = 0. This velocity field was advocated by Duchon and 
Robert |DRj as an example of 'relaxed solution' of the Euler equations. 
Notice the apparent separation of space variables: v± depends only on X\ 
(through £ = xit~ 2 / 3 ), v 3 depends only on x 3 while t> 2 = 0. Strictly speaking, 
this is not true, since formulae f)23|24p depend on the sign of x 3 — i^. However 
this is good enough to provide a very simple structure to the corresponding 
flow X(t,x) (defined by jTS): 

X(t,x) = (X 1 (t,Xi),x 2 ,X 3 (t,x 1 ,x 3 )), 

where the first component depends only on the first space variable (which is 
not the case of the last component). More precisely, by integration of (I23p . 
we get the following explicit formula: 

Xi(t,zi) = t 2/3 (2 v / xit- 2 / 3 - 1), 0<xi<t 2/3 , (25) 
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X,(t, Xl ) = t 2/3 (l - 2^- Xl t- 2 / 3 ), -t 2 / 3 < Xl <0, (26) 
Xi(t,xi) = xi, |xi|>t 2/3 . (27) 

Since v is divergence free and parallel to dD, X(t, ■) is a volume preserving 
map of D. As a consequence; X \ — > X 1 (t,x 1 ) must be a Lebesgue measure 
map of the interval [— L,L]. Indeed, for each continuous function f( X ), 

f( X )d X = / f(X(t, X ))d X = / f(Xi(t, xi), 22, A 3 (t, xi, X 3))d X id X 2d X s 

D JD JD 

and, in particular, when / = f( X ±): 



I f( Xl )d Xl = f f{X 1 (t, Xl ))d Xl . 

J-L J~L 



From the definition of V\, we also deduce 

d t + d(^-) + d lP = 0, (28) 

(in the sense of distribution, with no spurious singular measure), where 

P(t,Xi) = -^ f2 (t 4/3 -x 2 1 ) + . (29) 

Thus (t>i, t>2, t>3,p) solves the hydrostatic equations (in distribution form). 
Also notice that (1281) just means: 

a&x 1 (t,x 1 )) = -(d lP ){t,x 1 (t,xi)). 

Finally, by setting 

g t ( X ) = (X 1 (t, Xl ), X2 ,x 3 ), p t (x) =p(t, Xl ), 
we recover the solution discussed in Theorem 14.11 

Acknowledgments 

The author acknowledges the support of ANR contract OTARIE ANR-07- 
BLAN-0235. Part of his research was done during the a stay at UBC, Van- 
couver. (Thematic program and PDE summer school, workshop on Regu- 
larity problems in hydrodynamics, organized by V. Sverak and T.-P. Tsai, 
PIMS-UBC August 10-14, 2009.) 



12 



The minimizing geodesic problem in fluid mechanics 



References 

[AF1] L. Ambrosio, A. Figalli, Geodesies in the space of measure-preserving 
maps and plans, Arch. Rat. Mech. Anal. 194 (2009) 421-462. 

[AF2] L. Ambrosio, A. Figalli, On the regularity of the pressure field of Bre- 
mer' s weak solutions to incompressible Euler equations, Calc. Var. 
Partial Differential Equations 31 (2008) 497-509. 

[AK] V. I. Arnold, B. Khesin, Topological methods in Hydrodynamics, 
Springer Verlag, 1998. 

[Ar] V. I. Arnold, Mathematical methods of classical Mechanics, Graduate 
texts in mathematics, 60, Springer-Verlag, New-York, 1989 . 

[BFS] M. Bernot, A. Figalli, F. Santambrogio, Generalized solutions for the 
Euler equations in one and two dimensions, J. Math. Pures Appl. 91 
(2009) 137-155. 

[Brl] Y. Brenier, The least action principle and the related concept of gen- 
eralized flows for incompressible perfect fluids, J. of the AMS 2 (1989) 
225-255. 

[Br2] Y. Brenier, The dual least action principle for an ideal, incompressible 
fluid Arch. Rational Mech. Anal. 122 (1993) 323-351. 

[Br3] Y. Brenier, Minimal geodesies on groups of volume-preserving maps, 
Comm. Pure Appl. Math. 52 (1999) 411-452. 

[Br4] Y. Brenier, Homogeneous hydrostatic flows with convex velocity pro- 
files, Nonlinearity 12 (1999) 495-512. 

[Br5] Y. Brenier, Remarks on the derivation of the hydrostatic limit of the 
Euler equations Bull. Sci. Math. 127 (2003) 585-595. 

[Br6] Y. Brenier, Generalized solutions and hydrostatic approximation of the 
Euler equations, Physica D 237 (2008) 1982-1988. 

[BG] Y. Brenier, W. Gangbo, L p approximation of maps by diffeomor- 
phisms, Calc. Var. 16 (2003) 147-164. 

[DR] J. Duchon, R. Robert, Relaxation of the Euler equations and hydrody- 
namic instabilities, Quarterly Appl. Math. 50 (1992) 235-255. 



13 



The minimizing geodesic problem in fluid mechanics 



[EM] D. Ebin, J. Marsden, Groups of diffeomorphisms and the notion of an 
incompressible fluid, Ann. of Math. 92 (1970) 102-163. 

[Eu] L. Euler, Opera Omnia, Series Secunda, 12, 274-361. 

[Gr] E. Grenier, On the derivation of homogeneous hydrostatic equations, 
M2AN Math. Model. Numer. Anal. 33 (1999) 965-970. 

[Li] P. -L. Lions, Mathematical topics in fluid mechanics. Vol. 1. Incom- 
pressible models, Oxford Lecture Series in Mathematics and its Appli- 
cations, Oxford University Press, New York, 1996. 

[Ne] Y. Neretin, Categories of bistochastic measures and representations of 
some infinite- dimensional groups, Sb. 183 (1992), no. 2, 52-76. 

[Shi] A. Shnirelman, On the geometry of the group of diffeomorphisms and 
the dynamics of an ideal incompressible fluid, Math. Sbornik USSR 
56 (1987) 79-105. 

[Sh2] A. I. Shnirelman, Generalized fluid flows, their approximation and 
applications, Geom. Funct. Anal. 4 (1994) 586-620. 



14 



